Quantum Field Nature of Relaxation by Strong Coupling in Joint Electromagnetic 
Field - Dipole Moments - Phonon Field System 
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QED-model for dynamics of spectroscopic transitions in ID multiqubit exchange coupled system, 
presented in [l4| is generalized by taking into account more carefully quantum nature of relaxation 
processes. The formation of coherent system of the resonance phonons is predicted and experimen- 
tally confirmed. It leads to oscillation character of quantum relaxation, that was experimentally 
identified for the first time. 
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I. INTRODUCTION AND BACKGROUND 



Quantum electrodynamics (QED) takes on more and 
more significance for its practical application and it, in 
fact, becomes to be working instrument in spectroscopy 
studies and industrial spectroscopy control. Quantum 
dynamics of two level systems (qubits), coupled to a sin- 
gle mode of an electromagnetic cavity are of consider- 
able interest in connection with quickly developing new 
quantum physics branches like to cavity quantum elec- 
trodynamics [l| or quantum computing [2J. Spin-based 
solid state quantum bits are known to have long coher- 
ence times, while also offering the promise of scalabil- 
ity, and are natural building blocks for quantum com- 
putation. Phosphorus donor nuclei in silicon have been 
known since the 1950s to have some of the best spin co- 
herence properties in solids. The spin coherence time T2 
measured by Hahn spin echo method for donor electron 
spins in bulk Si:P has been reported to be equal ~ 60ms 
Q. This is the longest coherence time measured in elec- 
tron spin qubits, and greatly exceeds the value reported, 
for instance, in GaAs quantum dots, which, measured 
also by Hahn spin echo method, is equal to ~ 100/^s [|J. 
However, fabrication of ordered and gated donor arrays 
and coherent control over donor electrons has turned out 
to be extremely difficult. There is now much interest 
in fabricating spin-based devices in diamond, with po- 
tential applications in quantum communication, quan- 
tum computation, and magnetometry. Nitrogen- vacancy 
(NV) centers located deep within a diamond lattice ap- 
pear promising solid-state spin qubits since they com- 
bine optical initialization and readout capabilities owing 
to long electron spin coherence life times approaching 
lms at room temperature and the ability to control 
coupling to individual nuclear spins. In all of aforesaid 
applications it is necessary or at least advantageous to 
couple NV centers to optical structures like to waveg- 
uides and microresonators, to enable communication be- 
tween distant qubits or to allow efficient extraction of 
emitted photons. Therefore, a reliable method is needed 
to create NV centers with good spectral properties in 



close proximity (< lOOnm) to a diamond surface. In 
addition, the charge state of the NV center must be con- 
trolled. Given tasks are under study at present. We will 
show in present work another way to obtain long-lived 
coherent states with similar field of practical application. 
The prediction will be based on QED-theory of spectro- 
scopic transitions. The simplest models which capture 
the salient features of the relevant physics in given field 
are the Jaynes-Cummings model (JCM) Q for the one 
qubit case and its generalization for multiqubit systems 
by Tavis and Cummings Q • Tavis-Cummings model was 
generalized in 0, by taking into account the ID-coupling 
between qubits. Recently QED-model for one chain cou- 
pled qubit system was generalized for quasionedimen- 
sional axially symmetric multichain coupled qubit sys- 
tem Q. It is substantial, that in the model, proposed 
in H the interaction of quantized EM-field with multi- 
chain qubit system is considered by taking into account 
both the intrachain and interchain qubit coupling with- 
out restriction on their values. It follows from theoret- 
ical results injsj], Q and from their experimental con- 
firmation in [lj], that by strong interaction of EM-field 
with matter the correct description of spectroscopic tran- 
sitions including stationary spectroscopy is achieved the 
only in the frame of QED consideration. It concerns both 
optical and radio spectroscopies, that means, that QED 
consideration has to be also undertaken by electron spin 
resonance (ESR) studies in the case of strong interaction 
of EM-field with spin systems. It is reasonable to sug- 
gest, that analogous conclusion can be drawn for the case 
of strong interaction of phonon field with spin system or 
electron system. In other words it seems to be reasonable 
the idea, that relaxation of paramagnetic (or optical) cen- 
ters in the case of strong spin-phonon (electron-phonon) 
interaction can be described correctly the only in the 
frames of quantum field theory. 

The aim of the work presented is to show theoreti- 
cally and experimentally, that in the case aforesaid the 
relaxation process of excited matter electronic or matter 
spin subsystem by resonance interaction with EM-field, 
in particular by ESR studies, can really display its quan- 
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turn nature and to play the essential role in dynamics of 
spectroscopic transitions, including radiospectroscopy. 



is vector-operator of spectroscopic transitions for Ith 
chain unit, I = 2,7V- 1 [lj 
the operators 



Its components, that is, 



II. RESULTS AND DISCUSSION 

A. QED Equations for Spectroscopic Transitions 
by Strong Electron-Phonon Coupling 

Recently in the work [l4| the system of QED difference- 
differencial equations for dynamics of spectroscopic tran- 
sitions for both radio- and optical spectroscopy for the 
model, representing itself the ID-chain of N two-level 
equivalent elements interacting with quantized EM-field 
(that is N qubits) , coupled by exchange interaction (or its 
optical analogue for the optical transitions) have recently 
been derived. Naturally the equations are true for any 
3D system of paramagnetic centers (PC) or optical cen- 
ters by the absence of exchange interaction. In given case 
the model presented differs from Tavis-Cummings model 
[?J by inclusion into consideration of quantized phonon 
system, describing the relaxation processes from quan- 
tum fied theory position. Seven equations for the seven 
operator variables, describing joint system {field + mat- 
ter} can be presented in matrix form by three matrix 
equations. They are the following 
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ai m ee \j v )(m v \ (5) 

are set up in correspondence to the states \j v )^(m v \, 
where v = 1,N, j = a, /3, m = a, ft. For instance, the 
relationships for commutation rules are 
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where vector operators Gi-ij+i, I = 2, N — 1, are given 
by the expressions 
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In relations (0) J is the exchange interaction constant in 
the case of magnetic resonance transitions or its optical 
analogue in the case of optical transitions, the function 



V lk - in (TO]) is 
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where pj m is matrix element of operator of magnetic 
(electric) dipole moment Pi of l-th chain unit between 
the states \ji) and \mi) with j £ {a,/3}, m € {a,/?}, 
j =/= to, is unit polarization vector, ep is unit vec- 
tor along indirection, £^ is the quantity, which has the 

dimension of magnetic (electric) field strength, k is quan- 
tized EM-field wave vector, the components of which get 
a discrete set of values, ui^ is the frequency, corresponding 

to fcth mode of EM-field, ai and ag are EM-field creation 
and annihilation operators correspondingly. In the sug- 
gestion, that the contribution of spontaneous emission is 
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relatively small, we will have pj 7 



Pi 



j G {a, /3}, m £ {a,/3}, j 7^ m. Further matrix ||-R^ || is 



pi, where 

(A) 1 



1^1 



1 



ql I 



(12) 



Here B^' is 
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bi (6g-) is the creation (annihilation) operator of the 
phonon with impulse q and with energy fno^, is 
electron-phonon coupling constant. In equations ([2]) and 
([3]) 1 1 (Tp| j is Pauli z-matrix, in equation (fTJ) is diagonal 
matrix, numerical values of its elements are dependent on 
the basis choice. It is at appropriate basis 



themselves vector-operator difference-differential gener- 
alization of the system, which belongs to well known fam- 
ily of equation systems - Volterra model systems, widely 
used in biological tasks of population dynamics studies, 
which in its turn is generalization of Verhulst equation. 
In other words, it was predicted, for instance, that by 
some parameters in two-sybsystem Volterra model the 
stochastic component in solution will be appeared. Given 
prediction has experimental confirmation by the study 
of optical properties in carbynes [l6| indicating, that in 
given material strong electron-photon interaction is re- 
alized, which allows to explain the posibility to observe 
the stationary IR-rcflection or absorbtion spectra both in 
usual and in stochastic regime. 

The terms like to right hand side terms in ^ were 
used in so called "spin-boson" Hamiltonian [l 71 ] and in so 
called " independent boson model" 18]. Given models 
were used to study phonon effects in a single quantum 
dot within a microcavity F[l, Hj, [Hj], 0]> Si- So, 
it has been shown in 1231 , [23| , that the presence of the 
term in Hamiltonian 1 1411 
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Right hand side expression in ([7]) is vector product 
of vector operators. It can be calculated by using of 
known expression (II 5[) with additional coefficient ^ the 
only, which is appeared, since 
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the products of two components of two vector operators 
are replaced by anticommutators of corresponding com- 
ponents. Given detail is mapped by symbol <g) in (JT]) and 
by symbol ' in determinant (|15|) . 

It follows from comparison with semiclassical Landau- 
Lifshitz (L-L) equation for dynamics of spectroscopic 
transitions in a chain of exchange coupled centers [15], 
[l4j . that the equation, which is given by (fTJ) is its 
QED-generalization. In comparison with semiclassical 
description, where the description of dynamics of spec- 
troscopic transitions is exhausted by one vector equa- 
tion (L-L equation or L-L based equation), in the case 
of completely quantum consideration L-L type equation 
describes the only one subsystem of three-part-system, 
which consist of EM-ficld, dipole moments' (magnetic or 
electric) matter subsystem and phonon subsystem. It 
was concluded in [lj], that the presence of additional 
equations for description of transition dynamics by QED 
model in comparison with semiclassical model leads to a 
number of new effects, which can be predicted the only 
by QED consideration of resonance transition phenom- 
ena. One of new effect was described in starting 
the only from the mathematical structure of the equa- 
tions. It was argued, that the equations (flj, (0) represent 
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which coincides with corresponding term in Hamiltonian 
in [22j], [23| at N = 1 [contribution of given term to 

N 

the equations for spectroscopic transitions is &i^q, 

1=1 

see equation ([3]), (note that the equations for spec- 
troscopic transitions were not derived in above cited 

works [H, HI, EH. El, HI)] leads the onl y t0 ex- 
ponential decrease of the magnitude of quantum Rabi 
oscillations with increase of electron-phonon coupling 
strength and even to their supression at relatively strong 
electron-phonon coupling. However by strong electron- 
photon coupling amd strong electron-phonon coupling 
quite other picture of quantum relaxation processes be- 
comes to be possible. Really, if to define the wave func- 
tion of the chain system, interactig with quantized EM- 
field and with phonon field, to be vector of the state in 
Hilbert space over quaternion ring, that is quaternion 
function of quaternion argument, then like to |1J] can be 
shown, that the equations (flj to are Lorentz invari- 
ant and the transfer to obscrvablcs can be realized. In 
particular, taking into account, that quaternion vector of 
the state is proportional to spin, the Hamiltonian, given 
by (JTHJ) describes in fact the interaction of phonon field 
with z- component S z of the spin of matter subsystem. 
It seems to be reasonable to take into consideration the 
interaction of phonon field with S + - and S~ components 
of the spin of matter subsystem. Therefore we come in a 
natural way to the following Hamiltonian 
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where H c is chain Hamiltonian by the absence of the in- 
teraction with EM-field, / H F is field Hamiltonian, H CF 
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is Hamiltonian, describing the interaction between quan- 
tized EM-ficld and atomic chain. Hamiltonian T-L c is 



n c = n Q + n J , 



(18) 



where T-L° is chain Hamiltonian in the absence of the inter- 
action between structural elementary units of the chain. 
H is given by the expression 



n° = E mv \m v ) (m v \ 



(19) 



where at is EM-field creation operator, a r is EM-field 

k k 

annihilation operator, superscript * in Q*j% means com- 
plex conjugation. Field Hamiltonians have usual form 



for EM-held and 



(24) 



(25) 



Here m = a, /3, E mv are eigenvalues of "H , which corre- 
spond to the states \m v ) of vth chain unit. Hamiltonian 
n J is 

(20) 



for phonon held. Hamiltonian V. is 

njCPh _ njCPh , njCPh 
TL — t~L z t" rt^ , 



N 1 . 

H J = y^[J B (o-^o- ra+1 + <7„ (J^+i + Q^n^n+i) + H.c.]. where H^ ph is determined by the expression 



(26) 



It is suggested in the model, that \a n ) and |/3„) are eigen- 
states, producing the full set for each of N elements. It 
is evident, that given assumption can be realized strictly 
the only by the absence of the interaction between the 
elements. At the same time proposed model will rather 
well describe the real case, if the interaction energy of ad- 
jacent elements is much less of the energy of the splitting 
hojQ = £p — E a between the energy levels, corresponding 
to the states \a n ) and |/3 n ). The case considered includes 
in fact all known experimental situations. It is clear, 
that Hamiltonian 1-L of interaction of quantized EM- 
held with atomic chains can also be represented in the set 
of variables, which includes the components of spectro- 
scopic transition vector operator a v . Really, suggesting 
dipole approximation to be true and polarization of held 
components to be fixed, we have 
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j=l l^tm m J, (21) 
e -i"Rt+ifc + H.C.], 

where p 1 ™ 1 is matrix element of operator of magnetic 

(electric) dipole moment Pj of j-th chain unit between 
the states and \rrij) with lj — oej,0j,mj = otj,0j, e? 

is unit polarization vector, ep is unit vector along Pj- 
direction, (Sg is the quantity, which has the dimension of 

magnetic (electric) field strength, k is wave vector, is 
field annihilation operator. In the suggestion, that the 
contribution of spontaneous emission is relatively small, 
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Then the expression (|2ip can be rewritten in the form 
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Hamiltonian H± ph can be represented in the following 
form 



/V 



H 



CPh 



= EE^ 



(28) 

Here A^ and A- are electron-phonon coupling constants, 
which characterisire correspondingly the interaction with 
z- component Sj and with S + - and S~ components of 
the spin of jth chain unit. It seems to be understand- 
able, that they can be different in general case. More- 
over, in order to take into account the interaction with 
both equilibrium and nonequilibrium phonons both the 
electron-phonon coupling constants have to be complex 
numbers, that takes proper account by expressions (|2T[> . 

(EH). 

It can be shown, that the equations of the motion for 
spectroscopic transition operators ct; , for quantized EM- 
field operators ag, at and for phonon field operators b^, 

bt are the following. Instead equation (JTJ) we have 
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The equation © remains without changes. The equation 
© is 
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Thus, QED-model for dynamics of spectroscopic tran- 
sitions in ID multiqubit exchange coupled system, pre- 
sented in [141 ] is generalized by taking into account more 
carefully quantum nature of relaxation processes. New 
quantum phenomenon can be predicted. The prediction 
results from the structure of the equations derived and 
it means, that by resonance the coherent system of the 
resonance phonons, that is, the phonons with the en- 
ergy, egualed to resonance photon energy can be formed, 
which can lead to appearance along with Rabi oscillations 
determined by spin (electron)-photon coupling with the 
frequency £l RF of Rabi oscillations determined by spin 
(electron)-phonon coupling with the frequency fl RPh . In 
other words QED-model predicts the oscillation charac- 
ter of quantum relaxation, that is quite different charac- 
ter in comparison with phcnomcnological Bloch model. 
Moreover if < g the second Rabi oscillation process 
will be observed by stationary state of two subsystems 
{ EM-fied + magnetic (electric) dipoles }, that it will 
be registered in quadrature with the first Rabi oscilla- 
tion process. It can be experimentally detected even by 
stationary spectroscopy methods. 



B. Experiment 

ESR studies in anthracite samples have been carried 
out on stationary ESR spectrometer " Radiopan" at room 



temperature. The 100 kHz high frequency modulation 
of static magnetic field Hq was used. The static mag- 
netic field Hq was sweeped in two directions from lesser 
value to bigger value {^^- > 0) and from bigger value 
to lesser value (^f- < 0). Two regime of registration 
were used - with automatic microwave frequency adjust- 
ment and without automatic microwave frequency ad- 
justment. The spectra obtained are presented in Fig- 
ures 1, 3-9. The most interesting, that the spectra can 
be registered by using of 100 kHz high frequency (HF) 
modulation canal at modulation amplitude H m = 0, Fig- 
ure 1. It is seen from Figure 1, that the spectrum has 
very unusual for stationary ESR-spectroscopy form, it is 
similar partly to the spectra, registered by nonstation- 
ary transient ESR-spectroscopy. At the same time the 
functional dependence is another. The spectral distri- 
bution of nutation signals in transient ESR-spectroscopy 
is described usually by Bessel functions. The spectrum, 
presented in Figure 1 , consist of two group of very narrow 
ESR lines. For instance, the linewidth AH of the cen- 
tral line with maximal amplitude in left group is equal to 
0.01 (±0.004) G. It allows to determine g- value (that is 
g-factor in Zeeman term of spin-Hamiltonian) very pre- 
cisely. It is equal to 2.002735 (±0.000001). The linewidth 
of the other lines is practically the same. The lines pre- 
sented in Figure 1 seem to be the most narrow among all 
the lines, registered in ESR-spectroscopy at all. The po- 
sition of the central line with maximal amplitude in right 
group is given by g-value 2.00255 (±0.00001). The struc- 
ture of the spectrum is very similar to theoretical time 
evolution for a system with one qubit in JCM, Figure 2. 
Really, the spectra in both the groups can be described 
by almost sinusoidal dependence with Gaussian envelope, 
Figure 3. The qualitative difference between time evolu- 
tion given by JCM and spectral distribution presented in 
Figure 1 is that, the left group of lines in JCM has cosi- 
nusoidal dependence, Figure 2. The possibility to explain 
the spectrum observed by usual amplitude modulation is 
excluded, since the envelope by amplitude modulation is 
described by sinusoidal dependence, we have the enve- 
lope, described by Gauss distribution with slightly dif- 
ferent parameters for left and right sides. So the fit of 
experimental data, Figure 3, for left hand side group of 
lines gives the following expression 
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y = y° + w\ln eM ~ 



(H - H c ) 
W 



h 



(35) 



where y = 52.4, H c = 3322.95G, W = 0.36G, A = 32G. 

Given result and comparison with Figure 2 allow to 
conclude that the spectrum in Figure 1 represents itself 
the quantum Rabi oscillation picture. Really, let us touch 
on JCM in some details. Central role in JCM plays the 
sum with infinite summation limit, which represents on 
a scale [-1,1] the degree of excitation of two-level system 
resulting of interaction between the atom dipole or spin 
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Figure 1: ESR spectrun of anthracite sample, registered by Figure 3: The fit of left hand side envelope of left oscillation 
100 kHz high frequency modulation of static magnetic field group in the spectrum, presented in Figure 1 by Gaussian 
with amplitude H m = 
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Figure 2: Theoretical time evolution for a system with one 
qubit in JCM, initial state \ip)(g> \a) is direct product of atomic 
ground state and coherent field state with n — 50 
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Figure 4: The fit of the splitting between the lines in left 
oscillation group in the spectrum, presented in Figure 1 by 
linear dependence 



and single mode quantized EM-field. It is 

oo I 1 2n 

(a z (t)) =-exp[-|a| 2 ]^^Lcos2 ffv ^, (36) 

n=l U ' 

where \a) is fully coherent state of the field, taking place 
at t = 0, at that \a\ 2 = n, that is, it is the average 
number of photons in the field, g is coupling constant 
between field and atom (spin). The sum in Q36p cannot 
be expressed exactly in analytical form. For very short 
times and very large n behavior of (a z (t)) is determined 
by cos2gy/fit. Cummings [24| has shown, that by reso- 
nance and for intermediate time t values the cosine Rabi 
oscillation damp quickly (so called collapse takes place). 



Given damping can be described by Gaussian envelope 

eM-\(9tf]- (37) 

It is substantial, that it not depends on field intensity 
unlike to semiclassical Rabi oscillation damping process 
and it is determined entirely the o nly by coupling con- 
stant g. The authors of the work [25| have found, that 
JCM contains so called revival process with revival time 
Tr, given by the expression 

T fl = ^VA 2 +Vn, (38) 

g 2 

where A is deviation of field mode frequency from res- 
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Figure 5: ESR spectrun of anthracite sample, registered by 
100 kHz high frequency modulation of static magnetic field 
with amplitude H m = 0.01G and with automatic microwave 
frequency adjustment, 2ga < q 



Figure 7: ESR spectrun of anthracite sample, registered by 
100 kHz high frequency modulation of static magnetic field 
with amplitude Hm = 0.01G and with automatic microwave 
frequency adjustment, > 
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Magmetic Field (G) 

Figure 6: Central part of ESR spectrun of anthracite sample, 
registered by 100 kHz high frequency modulation of static 
magnetic field with amplitude H m — 0.01G and with auto- 
matic microwave frequency adjustment, < 



onance value. Revival process takes place at all time 
values, satisfying the relation t = kTn, k £ N. It is 
seen from (|38|) that revival time depends on n and it is 
proportional to field amplitude at A = like to depen- 
dence of Rabi frequency by semiclassical consideration. 
Let us illustrate the dynamics of spectroscopic transi- 
tions in JCM. The initial state for single qubit system 
can be defined by direct product of the two level matter 
(atomic, spin and so on) subsystem state \tp) and coher- 
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Figure 8: Central part of ESR spectrun of anthracite sample, 
registered by 100 kHz high frequency modulation of static 
magnetic field with amplitude H m — 0.01G and with auto- 
matic microwave frequency adjustment, > 

ent EM-field state It is the following 

|*(0)) = |$ ® \a), (39) 

where is 

|V) = ci|Vi)+C2hfe) (40) 

and | a) is 

I 1 2 00 I In 

\ a )=expf^]J2^ r \n),a = y/ne i *. (41) 

n— 1 
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Figure 9: ESR spectrun of anthracite sample, registered by 
100 kHz high frequency modulation of static magnetic field 
with amplitude H m — 0.01G and without automatic mi- 
crowave frequency adjustment, > 

The Rabi oscillations of the probability, that the qubit is 
in the initial state in the first damping stage (collapse), 

on a time scale of t c ~ — and then revive at Tr ~ 2ZLy5I 
g x g 

are illustrated in Figure 2 for c\ — 1, C2 = by plotting 

oo 

j^i&iMnm 2 , (42) 

n=l 

where (ipi,n\ = (ipi \ ® (n| is the ground state for the 
qubit with n photons in the cavity. Our measurements 
of the dependence of the ESR spectral characteristics on 
the microwave power (photon number) is also showed 
that in correspondence with JCM the positions of the 
lines of left hand side group are not dependent on pho- 
ton number, at the same time the positions of the lines 
of right hand side group are dependent on photon num- 
ber. The aforesaid difference between our results and 
JCM can be easily explained if we determine the Rabi 
frequency value. It is equaled to 28 kHz. Given value is 
substantially lesser in the same microwave power range, 
than the Rabi frequency values in excited by microwave 
power systems. So, for instance, in Si the Rabi frequency 
changes from ss 1 x 10 6 rads~ 1 to ps 4.5 x 10 6 rads~ 1 
for the Si-P3 centers (26|. It means that we really have 
observed the quantum relaxation oscillations, predicted 
theoretically in Section 2A. The prediction of sinusoidal 
character (instead cosinusoidal in JCM) is also confirmed. 
The comparison with the model for dynamic of spectro- 
scopic transitions, given in Section 2 A is correct, since at 
modulation amplitude H m = the absorption cannot be 
registered. At the same time the emission can be regis- 
tered, since by detection at 100 kHz the quantity which 
changes essentially slowly is registered like to usual ESR 
detection. Then, if to omit the part from Hamiltonian, 



corresponding to EM-held and z-relaxation term and to 
set N equaled to 1 , we will have in rotating wave approx- 
imation JCM Hamiltonian for phonon subsystem. The 
linear dependence of splitting, presented in Figure 4, can 
be explained by the dependence of Vt RPh on the static 
H-field, since by the change of static H-field the distance 
between the energy levels in qubits is also changes. Ex- 
panding the dependence £l RPh (H) in Taylor series and 
restricting by linear term (it is correct, since H-change is 
small), we obtain the agreement with experiment. 

To confirm the conclusion on coherent emission char- 
acter (which is like to maser emission character) of the 
spectrum, observed by H m = and to establish the struc- 
ture of the centers, which can emit stationary, we have 
undertaken the study of ESR-responce at various values 
of HF-modulation amplitude. The spectra, registered by 
100 kHz HF-modulation of static magnetic field with am- 
plitude H m — 0.01G are presented in Figures 5 to 9. 
The lines have also unusual shape for absorption deriva- 
tive. The spectra are also very different by the change of 
sweep direction, compare Figures 5 and 7, 6 and 8, which 
all were registered with automatic microwave frequency 
adjustment. Given difference disappears by registration 
without automatic microwave frequency adjustment and 
for both sweep directions we obtain the spectrum, pre- 
sented in Figure 9. It is substantial, that observed polar- 
ity of derivative of ESR responce by usual registration 
with automatic microwave frequency adjustment does 
not depend on sweep direction. At the same time, when 
ESR spectra by ^ > and by ^ < are absorbtion 
spectra the polarities of qiven spectra have to be different 
[271 ] . Consequently one of the spectrum is emission spec- 
trum. It is evident, taking into consideration the energy 
conservation law, that emission is observed by ^ < 0. 
Although the emission is coherent, it is not maser effect, 
since the amplification is absent (it seems to be in princi- 
ple impossible for two level system). Then the hysteresis 
effect (compare Figure 5 and 7), which is observed by 
registration with automatic microwave frequency adjust- 
ment can be explained in a natural way like to Stoks 
luminescence, where the emission spectrum is shifted to 
more low frequences in comparison with excitation spec- 
trum to be consequence of interaction with nonresonance 
phonons. The key to structure of ESR centers can be 
found by the analysis of the central part of ESR spectra 
by registration with automatic microwave frequency ad- 
justment, presented in Figures 6 and 8. It is seen from 
Figure 6, that Rabi oscillation process is also takes place, 
although spectral resolution is not very good, since the 
modulation amplitude is coinciding with linewidth. The 
most interesting is the shape of background line. In both 
the spectra the shape of background line has kink form 
and can be be approximated by tanh-function. It can 
mean, that the model of the ESR centers has to be sim- 
ilar to Su-Schrieffcr-Hccger (SSH) model of topological 
solitons in trans-polyacetylene [281 ] . The g- factor value, 
which is very near to g-factor value of neutral solitons in 
trans-polyacetylene testifies in favour of given proposal. 
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It can mean, that structural element of anthracite in- 
cludes carbon backbone of trans-polyacetylene. At the 
same time anthracite is structurally disordered material. 
It means, that the solitons has to be pinned and they 
have to give usual ESR spectra. To explain the appear- 
ance of kinks, at that it is substantionally, that they are 
characterized by different parameters for absorption and 
emission spectra, and the possibility to register Rabi os- 
cillations, we have to consider joint {EM-field + mag- 
netic dipole + resonance phonon field} system. EM-field 
and resonance phonon field produce radiation communi- 
cation in the direction of EM-field propagation and in 
time. It means that we have ID chain, at that the chains 
of two kinds are produced in time - with the sequence 
- absorption time segment - emission time segment and 
with reversal sequence. It is substantionally, that given 
time segments are different in their values. In fact the 
dimerisation take place like to trans-polyacetylene, how- 
ever instead space z-axis along t-axis. It is in agreement 
with known equality of rights of coordinates in Minkowski 
space. In other words the time has to be considered be- 
ing to be quantized in given case. The junction of given 
time ordered chains produces time topological violations 
of two kinds. The first topological violation can effec- 
tively absorb photons and emit resonance phonons, the 
second topological violation has reversal properties. So, 
we come to the model of instanton, which is like to some 
extent to SSH-soliton model. It is understandable, that 
the instantons of both kinds can free, that is without ad- 
ditional activation energy, move along t-axis like to SSH- 
solitons along space z-axis. It is the first experimental 
detection of instantons (to our knowledge). The known 
instantons in the literature are considered to be a special 
kind of vacuum oscillations in gluon field theory and also 
in gravitation field theory [29| . They have quite another 
origin. 

Therefore we have observed new kind of coherent 
states. They are modes of resonance phonon field. The 



coherence time value in the samples studied is evaluted 
to be equal at room temperature ~ 1 Orris, at that at rel- 
ativelly low microwave power and it will increase with 
power increase, achieving the maximal known values at 
present and even more. Substantial advantage for pos- 
sible applications is that, that the responce on external 
EM-field is very strong. At the same time the detection 
of single spin is rather technically difficult task. The H- 
field dependence allows to separate the phonon modes, 
the dependence on H-field direction seems to be also es- 
sential for practical applications. 



III. CONCLUSIONS 

QED-modcl for dynamics of spectroscopic transitions 
in ID multiqubit exchange coupled system, presented in 
[3] is generalized by taking into account more carefully 
quantum nature of relaxation processes. It is predicted, 
that by resonance the coherent system of the resonance 
phonons, that is, the phonons with the energy, equaled 
to resonance photon energy can be formed. It leads 
to appearance along with Rabi oscillations determined 
by spin (electron)-photon coupling with the frequency 
£l RF of Rabi oscillations determined by spin (electron)- 
phonon coupling with the frequency Sl R . Therefore 
QED-model predicts the oscillation character of quan- 
tum relaxation, that is quite different character in com- 
parison with phenomenological Bloch model. Moreover 
if Ajp < g the second Rabi oscillation process will be ob- 
served by stationary state of two subsystems {EM-fied 
+ magnetic (electric) dipoles}, that means, that it will 
be registered in quadrature with the first Rabi oscilla- 
tion process. It can be experimentally detected even by 
stationary spectroscopy methods. All the predictions are 
experimentally confirmed. The model of new kind in- 
stantons is proposed. 
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